Single nucleotide polymorphisms (SNPs) in the open reading frame (ORF) of prion protein gene (PRNP) in Nigerian livestock species.

BACKGROUND: Prion diseases, also known as transmissible spongiform encephalopathies (TSEs) remain one of the deleterious disorders, which have affected several animal species. Polymorphism of the prion protein (PRNP) gene majorly determines the susceptibility of animals to TSEs. However, only limited studies have examined the variation in PRNP gene in different Nigerian livestock species. Thus, this study aimed to identify the polymorphism of PRNP gene in Nigerian livestock species (including camel, dog, horse, goat, and sheep). We sequenced the open reading frame (ORF) of 65 camels, 31 village dogs and 12 horses from Nigeria and compared with PRNP sequences of 886 individuals retrieved from public databases. RESULTS: All the 994 individuals were assigned into 162 haplotypes. The sheep had the highest number of haplotypes (n = 54), and the camel had the lowest (n = 7). Phylogenetic tree further confirmed clustering of Nigerian individuals into their various species. We detected five non-synonymous SNPs of PRNP comprising of G9A, G10A, C11G, G12C, and T669C shared by all Nigerian livestock species and were in Hardy-Weinberg Equilibrium (HWE). The amino acid changes in these five non-synonymous SNP were all "benign" via Polyphen-2 program. Three SNPs G34C, T699C, and C738G occurred only in Nigerian dogs while C16G, G502A, G503A, and C681A in Nigerian horse. In addition, C50T was detected only in goats and sheep. CONCLUSION: Our study serves as the first to simultaneously investigate the polymorphism of PRNP gene in Nigerian livestock species and provides relevant information that could be adopted in programs targeted at breeding for prion diseases resistance.

Multiple Origins and Genomic Basis of Complex Traits in Sighthounds.

Sighthounds, a distinctive group of hounds comprising numerous breeds, have their origins rooted in ancient artificial selection of dogs. In this study, we performed genome sequencing for 123 sighthounds, including one breed from Africa, six breeds from Europe, two breeds from Russia, and four breeds and 12 village dogs from the Middle East. We gathered public genome data of five sighthounds and 98 other dogs as well as 31 gray wolves to pinpoint the origin and genes influencing the morphology of the sighthound genome. Population genomic analysis suggested that sighthounds originated from native dogs independently and were comprehensively admixed among breeds, supporting the multiple origins hypothesis of sighthounds. An additional 67 published ancient wolf genomes were added for gene flow detection. Results showed dramatic admixture of ancient wolves in African sighthounds, even more than with modern wolves. Whole-genome scan analysis identified 17 positively selected genes (PSGs) in the African population, 27 PSGs in the European population, and 54 PSGs in the Middle Eastern population. None of the PSGs overlapped in the three populations. Pooled PSGs of the three populations were significantly enriched in "regulation of release of sequestered calcium ion into cytosol" (gene ontology: 0051279), which is related to blood circulation and heart contraction. In addition, ESR1, JAK2, ADRB1, PRKCE, and CAMK2D were under positive selection in all three selected groups. This suggests that different PSGs in the same pathway contributed to the similar phenotype of sighthounds. We identified an ESR1 mutation (chr1: g.42,177,149 T > C) in the transcription factor (TF) binding site of Stat5a and a JAK2 mutation (chr1: g.93,277,007 T > A) in the TF binding site of Sox5. Functional experiments confirmed that the ESR1 and JAK2 mutation reduced their expression. Our results provide new insights into the domestication history and genomic basis of sighthounds.

Target capture sequencing for the first Nigerian genotype I ASFV genome.

African swine fever (ASF) is a contagious viral disease that affects domestic pigs and wild boars, causing significant economic losses globally. After the first Nigerian outbreak in 1997, there have been frequent reports of ASF in pig-producing regions in the country. To facilitate control, it is important to understand the genotype and phylogenetic relationship of ASF viruses (ASFVs). Recent genetic analysis of Nigerian ASFV isolates has revealed the presence of both genotypes I and II; this is based on analysis of a few selected genes. Phylogenetic analysis of ASFV whole genomes highlights virus origins and evolution in greater depth. However, there is currently no information on the ASFV genome from Nigerian isolates. Two ASFV-positive samples were detected during a random survey of 150 Nigerian indigenous pig samples collected in 2016. We assembled near-complete genomes of the two ASFV-positive samples using in-solution hybrid capture sequencing. The genome-wide phylogenetic tree assigned these two genomes into p72 genotype I, particularly close to the virulent Benin 97/1 strain. The two ASFVs share 99.94 and 99.92 % genomic sequence identity to Benin97/1. This provides insight into the origin and relationship of ASFV strains from Nigeria and Italy. The study reports for the first time the determination of near-complete genomes of ASFV using in-solution hybrid capture sequencing, which represents an important advance in understanding the global evolutionary landscape of ASFVs.

Population structure and evolutionary history of the greater cane rat (Thryonomys swinderianus) from the Guinean Forests of West Africa.

Grasscutter (Thryonomys swinderianus) is a large-body old world rodent found in sub-Saharan Africa. The body size and the unique taste of the meat of this major crop pest have made it a target of intense hunting and a potential consideration as a micro-livestock. However, there is insufficient knowledge on the genetic diversity of its populations across African Guinean forests. Herein, we investigated the genetic diversity, population structures and evolutionary history of seven Nigerian wild grasscutter populations together with individuals from Cameroon, Republic of Benin, and Ghana, using five mitochondrial fragments, including D-loop and cytochrome b (CYTB). D-loop haplotype diversity ranged from 0.571 (± 0.149) in Republic of Benin to 0.921 (± 0.013) in Ghana. Within Nigeria, the haplotype diversity ranged from 0.659 (± 0.059) in Cross River to 0.837 (± 0.075) in Ondo subpopulation. The fixation index (FST), haplotype frequency distribution and analysis of molecular variance revealed varying levels of population structures across populations. No significant signature of population contraction was detected in the grasscutter populations. Evolutionary analyses of CYTB suggests that South African population might have diverged from other populations about 6.1 (2.6-10.18, 95% CI) MYA. Taken together, this study reveals the population status and evolutionary history of grasscutter populations in the region.

Polymorphism of prion protein gene (PRNP) in Nigerian sheep.

Polymorphism of the prion protein gene (PRNP) gene determines an animal's susceptibility to scrapie. Three polymorphisms at codons 136, 154, and 171 have been linked to classical scrapie susceptibility, although many variants of PRNP have been reported. However, no study has investigated scrapie susceptibility in Nigerian sheep from the drier agro-climate zones. In this study, we aimed to identify PRNP polymorphism in nucleotide sequences of 126 Nigerian sheep by comparing them with public available studies on scrapie-affected sheep. Further, we deployed Polyphen-2, PROVEAN, and AMYCO analyses to determine the structure changes produced by the non-synonymous SNPs. Nineteen (19) SNPs were found in Nigerian sheep with 14 being non-synonymous. Interestingly, one novel SNP (T718C) was identified. There was a significant difference (P < 0.05) in the allele frequencies of PRNP codon 154 between sheep in Italy and Nigeria. Based on the prediction by Polyphen-2, R154H was probably damaging while H171Q was benign. Contrarily, all SNPs were neutral via PROVEAN analysis while two haplotypes (HYKK and HDKK) had similar amyloid propensity of PRNP with resistance haplotype in Nigerian sheep. Our study provides valuable information that could be possibly adopted in programs targeted at breeding for scrapie resistance in sheep from tropical regions.

Long amplicon HiFi sequencing for mitochondrial DNA genomes.

Long-read sequencing technology is a powerful approach with application in various genetic and genomic research. Herein, we developed the pipeline for long amplicon high-fidelity (HiFi) sequencing and then applied it for sequencing mitochondrial DNA (mtDNA) genomes from pools of 79 Tibetan Mastiffs. We amplified the mtDNA genome with long-range PCR using two pairs of primers. Two rounds of circular consensus sequencing (CCS) were conducted and their accuracy was evaluated. The results indicate that the second round of CCS can improve the accuracy of HiFi reads. In addition, the analysis of 79 high-quality mtDNA genomes shows the Tibetan Mastiffs from outside of the Tibetan Plateau experienced hybridization with other dogs. The high quality reads generator (HQGR) software is provided to facilitate data analyses, which is publicly accessible on GitHub (https://github.com/Caizf-script/HQGR). Our long amplicon HiFi sequencing pipeline can also be applied in various target enrichment strategies for small genomes and candidate genes.

Scrapie-associated polymorphisms of the prion protein gene (PRNP) in Nigerian native goats.

Scrapie is a fatal prion protein disease stiffly associated with single nucleotide polymorphism (SNPs) of the prion protein gene (PRNP). The prevalence of this deadly disease has been reported in small ruminants, including goats. The Nigerian goats are hardy, trypano-tolerant, and contribute to the protein intake of the increasing population. Although scrapie has been reported in Nigerian goats, there is no study on the polymorphism of the PRNP gene. Herein, we evaluated the genetic and allele distributions of PRNP polymorphism in 132 Nigerian goats and compared them with publicly available studies on scrapie-affected goats. We utilized Polyphen-2, PROVEAN and AMYCO programs to examine structural variations produced by the non-synonymous SNPs. Our study revealed 29 SNPs in Nigerian goats, of which 14 were non-synonymous, and 23 were novel. There were significant differences (P < 0.001) in the allele frequencies of PRNP codons 139, 146, 154 and 193 in Nigerian goats compared with scrapie-affected goats, except for Northern Italian goats at codon 154. Based on the prediction by Polyphen-2, R139S and N146S were 'benign', R154H was 'probably damaging', and T193I was 'possibly damaging'. In contrast, PROVEAN predicted 'neutral' for all non-synonymous SNPs, while AMYCO showed a similar amyloid propensity of PRNP for resistant haplotype and two haplotypes of Nigerian goats. Our study is the first to investigate the polymorphism of scrapie-related genes in Nigerian goats.

African Suid Genomes Provide Insights into the Local Adaptation to Diverse African Environments.

African wild suids consist of several endemic species that represent ancient members of the family Suidae and have colonized diverse habitats on the African continent. However, limited genomic resources for African wild suids hinder our understanding of their evolution and genetic diversity. In this study, we assembled high-quality genomes of a common warthog (Phacochoerus africanus), a red river hog (Potamochoerus porcus), as well as an East Asian Diannan small-ear pig (Sus scrofa). Phylogenetic analysis showed that common warthog and red river hog diverged from their common ancestor around the Miocene/Pliocene boundary, putatively predating their entry into Africa. We detected species-specific selective signals associated with sensory perception and interferon signaling pathways in common warthog and red river hog, respectively, which contributed to their local adaptation to savannah and tropical rainforest environments, respectively. The structural variation and evolving signals in genes involved in T-cell immunity, viral infection, and lymphoid development were identified in their ancestral lineage. Our results provide new insights into the evolutionary histories and divergent genetic adaptations of African suids.

Chromosome-level genome assembly of the Muscovy duck provides insight into fatty liver susceptibility.

The Muscovy duck (Cairina moschata) is an economically important poultry species, which is susceptible to fatty liver. Thus, the Muscovy duck may serve as an excellent candidate animal model of non-alcoholic fatty liver disease. However, the mechanisms underlying fatty liver development in this species are poorly understood. In this study, we report a chromosome-level genome assembly of the Muscovy duck, with a contig N50 of 11.8 Mb and scaffold N50 of 83.16 Mb. The susceptibility of Muscovy duck to fatty liver was mainly attributed to weak lipid catabolism capabilities (fatty acid ß-oxidation and lipolysis). Furthermore, conserved noncoding elements (CNEs) showing accelerated evolution contributed to fatty liver formation by down-regulating the expression of genes involved in hepatic lipid catabolism. We propose that the susceptibility of Muscovy duck to fatty liver is an evolutionary by-product. In conclusion, this study revealed the potential mechanisms underlying the susceptibility of Muscovy duck to fatty liver.

Genome-wide investigations reveal the population structure and selection signatures of Nigerian cattle adaptation in the sub-Saharan tropics.

BACKGROUND: Cattle are considered to be the most desirable livestock by small scale farmers. In Africa, although comprehensive genomic studies have been carried out on cattle, the genetic variations in indigenous cattle from Nigeria have not been fully explored. In this study, genome-wide analysis based on genotyping-by-sequencing (GBS) of 193 Nigerian cattle was used to reveal new insights on the history of West African cattle and their adaptation to the tropical African environment, particularly in sub-Saharan region.  RESULTS: The GBS data were evaluated against whole-genome sequencing (WGS) data and high rate of variant concordance between the two platforms was evident with high correlated genetic distance matrices genotyped by both methods suggestive of the reliability of GBS applicability in population genetics. The genetic structure of Nigerian cattle was observed to be homogenous and unique from other African cattle populations. Selection analysis for the genomic regions harboring imprints of adaptation revealed genes associated with immune responses, growth and reproduction, efficiency of feeds utilization, and heat tolerance. Our findings depict potential convergent adaptation between African cattle, dogs and humans with adaptive genes SPRY2 and ITGB1BP1 possibly involved in common physiological activities. CONCLUSION: The study presents unique genetic patterns of Nigerian cattle which provide new insights on the history of cattle in West Africa based on their population structure and the possibility of parallel adaptation between African cattle, dogs and humans in Africa which require further investigations.

Genetic diversity and population structure of muscovy duck (Cairina moschata) from Nigeria.

The domestic Muscovy duck (Cairina moschata) provide unique genetic resources patterned by both tropical environmental conditions and human activities, the evaluation of their genetic diversity and population structure will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the variation in mtDNA cytochrome b and nuclear DNA CYP2U1 sequences of 378 Nigerian Muscovy ducks (comprising of 287 de novo and 91 downloaded) plus 80 published sequences of Muscovy ducks from India. The results showed high haplotype diversity (0.800 ± 0.023) among Nigerian Muscovy duck populations with 91 distinct haplotypes for the nuclear DNA CYP2U1 gene but low (0.266 ± 0.033) for cytochrome b with 31 haplotypes. The median-joining networks of both markers grouped Nigerian Muscovy ducks into two; the first group consisting of only Nigerian Muscovy duck populations, and the second group Nigerian with Indian populations. Neutrality test results indicated that Nigerian populations experienced recent population expansion and/or genetic hitchhiking. A geographic signal was absent in line with previously studied poultry species in Nigeria. The most prominent haplotype dominated across all regions in Nigeria, which may be due to extensive genetic intermixing except for the Indian population (F ST = 0.02550, P = 0.01075). This indicated low genetic differentiation between and within Nigerian Muscovy duck as revealed by the suitability of the nuclear DNA CYP2U1 gene.

Genetic Diversity Relationship in Azakheli Buffalo Inferred from mtDNA and MC1R Sequences Comparison.

Azakheli is relatively smaller riverine breed with a very peculiar characteristics kept under unique traditional husbandry practices in comparison with rest of the Pakistani buffalo breeds; however, milk production is comparable. The present study was conducted to evaluate the genetic diversity of Azakheli breed. A total of sixty-six blood samples were collected for the amplification of mtDNA D-loop region and MC1R gene sequencing analysis. Median-joining network analysis of 191 mtDNA D-loop sequences of Azakheli and eight Indian riverine buffalo breeds clustered into three clades. Ancient Azakheli Region 1 clade was the oldest with the highest mutation steps and was present close to the root of UPGMA phylogenetic tree. There was 5 mutated lines distance between Pakistan buffalo and Indian riverine buffaloes. The populations of neighboring countries did not share any haplotypes with Azakheli buffalo of Pakistan. Possibly, residing for so long in the cold atmosphere and high elevation regions caused the mutation in mtDNA D-loop, though these conditions did not affect the overall performance of Azakheli as milch buffalo breed of Pakistan. MC1R analyses showed high mutations in Azakheli of Albino phenotype and all the black phenotype individuals of Azakheli buffalo share haplotypes with dominant Chinese and Indian black phenotypes buffaloes in MC1R median-joining network, indicating the reason of black coat color is due to MC1R gene. The haplotype diversity and nucleotide diversity was (H. 0.923, Pi: 0.00895) in Azakheli. Current results illustrated Asian ancestry for Azakheli buffalo, and mtDNA and MC1R analyses provided further evidence. Additional genetic analyses and archeological studies may provide further insight into the domestication period and history of Azakheli buffalo breed. The further studies are required on different coat colors with different genes on Azakheli buffalo to understand the phenotype variation.

Single-cell RNA Sequencing Reveals Thoracolumbar Vertebra Heterogeneity and Rib-genesis in Pigs.

Development of thoracolumbar vertebra (TLV) and rib primordium (RP) is a common evolutionary feature across vertebrates, although whole-organism analysis of the expression dynamics of TLV- and RP-related genes has been lacking. Here, we investigated the single-cell transcriptome landscape of thoracic vertebra (TV), lumbar vertebra (LV), and RP cells from a pig embryo at 27 days post-fertilization (dpf) and identified six cell types with distinct gene expression signatures. In-depth dissection of the gene expression dynamics and RNA velocity revealed a coupled process of osteogenesis and angiogenesis during TLV and RP development. Further analysis of cell type-specific and strand-specific expression uncovered the extremely high level of HOXA10 3'-UTR sequence specific to osteoblasts of LV cells, which may function as anti-HOXA10-antisense by counteracting the HOXA10-antisense effect to determine TLV transition. Thus, this work provides a valuable resource for understanding embryonic osteogenesis and angiogenesis underlying vertebrate TLV and RP development at the cell type-specific resolution, which serves as a comprehensive view on the transcriptional profile of animal embryo development.

Genomes reveal selective sweeps in kiang and donkey for high-altitude adaptation.

Over the last several hundred years, donkeys have adapted to high-altitude conditions on the Tibetan Plateau. Interestingly, the kiang, a closely related equid species, also inhabits this region. Previous reports have demonstrated the importance of specific genes and adaptive introgression in divergent lineages for adaptation to hypoxic conditions on the Tibetan Plateau. Here, we assessed whether donkeys and kiangs adapted to the Tibetan Plateau via the same or different biological pathways and whether adaptive introgression has occurred. We assembled a de novo genome from a kiang individual and analyzed the genomes of five kiangs and 93 donkeys (including 24 from the Tibetan Plateau). Our analyses suggested the existence of a strong hard selective sweep at the EPAS1 locus in kiangs. In Tibetan donkeys, however, another gene, i.e., EGLN1, was likely involved in their adaptation to high altitude. In addition, admixture analysis found no evidence for interspecific gene flow between kiangs and Tibetan donkeys. Our findings indicate that despite the short evolutionary time scale since the arrival of donkeys on the Tibetan Plateau, as well as the existence of a closely related species already adapted to hypoxia, Tibetan donkeys did not acquire adaptation via admixture but instead evolved adaptations via a different biological pathway.

Genomic Analyses of Unveil Helmeted Guinea Fowl (Numida meleagris) Domestication in West Africa.

Domestication of the helmeted guinea fowl (HGF; Numida meleagris) in Africa remains elusive. Here we report a high-quality de novo genome assembly for domestic HGF generated by long- and short-reads sequencing together with optical and chromatin interaction mapping. Using this assembly as the reference, we performed population genomic analyses for newly sequenced whole-genomes for 129 birds from Africa, Asia, and Europe, including domestic animals (n = 89), wild progenitors (n = 34), and their closely related wild species (n = 6). Our results reveal domestication of HGF in West Africa around 1,300-5,500 years ago. Scanning for selective signals characterized the functional genes in behavior and locomotion changes involved in domestication of HGF. The pleiotropy and linkage in genes affecting plumage color and fertility were revealed in the recent breeding of Italian domestic HGF. In addition to presenting a missing piece to the jigsaw puzzle of domestication in poultry, our study provides valuable genetic resources for researchers and breeders to improve production in this species.

Genetic Architecture Underlying Nascent Speciation-The Evolution of Eurasian Pigs under Domestication.

Speciation is a process whereby the evolution of reproductive barriers leads to isolated species. Although many studies have addressed large-effect genetic footprints in the advanced stages of speciation, the genetics of reproductive isolation in nascent stage of speciation remains unclear. Here, we show that pig domestication offers an interesting model for studying the early stages of speciation in great details. Pig breeds have not evolved the large X-effect of hybrid incompatibility commonly observed between "good species." Instead, deleterious epistatic interactions among multiple autosomal loci are common. These weak Dobzhansky-Muller incompatibilities confer partial hybrid inviability with sex biases in crosses between European and East Asian domestic pigs. The genomic incompatibility is enriched in pathways for angiogenesis, androgen receptor signaling and immunity, with an observation of many highly differentiated cis-regulatory variants. Our study suggests that partial hybrid inviability caused by pervasive but weak interactions among autosomal loci may be a hallmark of nascent speciation in mammals.